ON THE PROOF OF UNIVERSALITY FOR ORTHOGONAL AND 
SYMPLECTIC ENSEMBLES IN RANDOM MATRIX THEORY 
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O . 

■ Abstract. We give a streamlined proof of a quantitative version of a result 

from IDG1I which is crucial for the proof of universality in the bulk IDG1I and 
also at the edge |DG2| for orthogonal and symplectic ensembles of random 

O , matrices. As a byproduct, this result gives asymptotic information on a certain 

■ ratio of the j3 = 1,2,4 partition functions for log gases. 
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m— 1 
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(3) Q(q)=I(q) + ± l . 

For n = 2m — 1, define the (m — 1) x (m — 1) matrix 



( 4 ) T [m-1] = j _ H) 2 Q[m-l]_g[m-l] = j _ ^[m-1] 



where 

Qg 1 - 11 = Q( n - 2z + 2j), = 2m(. " f \ 1 < i, j<m- 1. 

Here Q = for jfe < 0. 

In |DG11 Theorem 2.6], the authors prove the following result: for m > 2, 

(5) detT^-^^O 

(see Remark |2 after the proof of Theorem ^ below) . Note that in the notation of 

pTf] . ri™- 1 ! = r m _L 

In this paper we will give a streamlined proof of the following quantitative version 
of ©. 

Theorem 1. For m > 2, 

(6) detTl" 1 - 1 ! > 0.0865. 

i 



2 



COSTIN, DEIFT AND GIOEV 



Equation JSJ plays a crucial role in proving universality in the bulk |DGlj . and 



also at the edge |DG2j . for orthogonal ((3 — 1) and symplectic (/3 = 4) random 
matrix ensembles for a class of weights w(x) = e^ v ^ where V(x) is a polyno- 
mial V(x) — K 2m i 2m + • • • , «2m > 0- (Here m is the same integer as in |JS}, |JB}.) 
The situation is as follows. In DGll IUG*2j . and also in |DGKV| . the authors use 
the method of Widom [W|, which is based in turn on |T W| . together with the as- 
ymptotic analysis for orthogonal polynomials in DKMVZ . A new and challenging 
feature of the method in [W] , which does not arise in the proof of universality in 
the case /3 = 2, is the appearance of the inverse of a certain matrix G\\ of fixed 
size n = 2m — 1 (see |DG1I (1.37) and Theorem 2.3 et seq.]). In the scaling limit 
as N — > oo, the matrix C\\ converges to a matrix and 

(7) detCft = (detTl 1 ™- 1 !) 2 



(see discussion from (2.13) up to Theorem 2.4 in |DGlj \ Thus in order to control 
the scaling limit for (3 = 1 and 4, we need to show that det T^™ -1 ! =/= 0. 

It turns out that detTl" 1-1 ! is related to partition functions for finite log gases 
in an external field V at inverse temperatures (3 = 1,2,4 



(8) 



ry / •• / II \x i -x i fe-^ v ^dx 1 ^dx k 



l<i<j<k 



1 1 ■ ■ ■ I e -^i<i<.<* Io K |« 1 -x,|-eJU v('t) dxi ... dx k . 



Using standard formulae for such partition functions (see e.g. |AvMI (4.4), (4.17), 
(4.20)]), together with [UGH (2.18)], one finds (see [SB Remark 2.4], [HUT! Remark 
1.5]) that for ensembles of (even) size N 

Thus 

Formula 0, together with {JJ, raises the possibility of using the methods of sta- 
tistical mechanics to prove iJSJ, (JSJ. The estimates in [J] show that the partition 
functions Zv,p,k have, for certain constants ay.p, leading order asymptotics of the 
form e a v,/3 k2 ( 1 +°( 1 )) as fc _^ qq^ anc i moreover, their combined contributions to 
detCn cancel to this order. In order to achieve cancellation at subsequent orders, 
and so prove (0), ©, one needs higher order asymptotics for the Zv,/3,k's, but, un- 
fortunately such asymptotics are known only for f3 = 2 (see |EM| ) . Regarding ijjj}, 
we take the contrary point of view, i.e., (|10fl and @ provide new information on 
the asymptotics of partition functions for log gases at inverse temperatures (3—1 
and 4. 

Much of the analysis in |DG1| involves estimating Q(q) in two regions: 3 < q < 
y/m and \fm < q < Am — 5. In this note, using bounds on 

(a) W{x) ^ir^3i ds 

7r Jo sin s 

which arc uniform in q = 3, 5, • • • and in < x < 7r/2 (see Lemma 0] below), we are 
able to estimate Q(q) uniformly for q = 3, 5, • • • , 4m — 5 and so avoid many of the 
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technicalities in the proof in |DG1| of (0- Of course the function W(x) is familiar 
from the analysis of the Gibbs phenomenon in Fourier analysis. 

Remark 1. For m = 1, corresponding to the Gaussian orthogonal and symplectic 
ensembles, T^ m ~^ is not defined and no analog of (JHJ), © is needed (see |DGlj ). 



We use the following result. For a matrix X let r(X) = sup{|A| : A € spec AT} 
denote the spectral radius of X . As is well known, for any operator norm |j ■ || on 
{X}, 

(12) r(X) = lim WX^ 1 ^ = inf H-X^ll 1 ^'. 

j-»oo j>l 

Lemma 2. Assume K and K' are J -dimensional matrices with real entries such 
that \Kij\ < K-p 1 < i,j < J, and r(K') < 1. Then r(K) < 1 and 

(13) det(7 - K) > det(J - K') > 0. 
Proof. The following is true: if r(X) < 1, then 

(14) det(7-X) = e -^™^ trix ' ) . 

This result is usually stated in the form that (|14fl holds if ||X|| < 1 (see e.g. |KeSi| ) . 
To obtain ljTi|) for r(X) < 1 from the case ||X|| < 1 simply apply (jT4*)) to [iX for 
H small and observe that for any fixed p satisfying r(X) < p < 1, < p l for I 

sufficiently large: then (|14fl follows for r(X) < 1 by analytic continuation fi —> 1. 

Equip M. J with the /oo-norm || • ||oo (any Z p -norm, 1 < p < oo would do) and for 
a matrix X mapping M. J — > R J denote the associated operator norm by \\X\\. For 
(j> = {4>j} G M. J we denote the vector with coordinates {|0j|} by \<j>\. We claim that 
r(K) < r(K'). Indeed, for S R J , \{K l ^)j\ < ((Jf')'l0l)j and so 

H^Vlloo < II {K') l \4>\ Woo < \\{K') l \\ || |0| |U = \\{K') l \\ ii^iu. 

Thus \\K l \\ < \\(K') l \\ and so r(K) < r(K') < 1 by (H2J). It follows that JUJ is valid 
for K and K' . But clearly | tr(K l )\ < tr((K') 1 ) and lO is now immediate. □ 

The function h(x) in has the following properties (see |DG11 Proposition 
6.2]): for < x < 1 

(i) h solves the differential equation 

x{x 2 - l)ti + (2m - 1 - 2(m - l)x 2 )h = 4m 

(15) (U) 4m = /i(0) < h(x) < h(l) = 4m 

2m — 1 

. . j / \ 4ml 2 " 1 - 1 f 1 t~ 2m , 
{Hi) h(x) = — ; / dt. 

Property (i) reflects the fact that h is a hypergeometric function, 

, , 4m , , 

K x ) = 71 72-Fi(l, -m + 1, -m + 3/2; x 2 ) 

2m — 1 

(see |DG1I (6.11)]) and (iii) follows by integrating (i). Property (ii) follows from (i) 
and CJ. 

Set 

(16) u(x ) Stt ( a;m ) = _l__i^L + i L. 
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Note that the function u(x) is closely related to the function y m which plays a 



prominent role in |DGlj : we have 



. . \/l — x 2 .1 

u(x) = y m (axcsmx) H , < x < 1. 

m 2m 

Also note that using the elementary identities for q = 3, 5, • • • , 

2 r /2 , /vr\ 2 r /2 sings , 

— / smgs sin s ds — 0, W I — I = — / — : as = 1 

Ti" Jo ^ 2 / 7T J sin s 

we have from J5J) , (|TB|) 

^ ri \ 4 r /2sin ^ / \ j 1 

(17) 7(g) = — / — u(coss)as— - — . 

n J sms 2m 

The main technical result in our proof of Theorem ^ is the following. 

Lemma 3. The function u(x) — u(x; m), m > 2, has the following properties. 

(i) u(x) is unimodal for x £ [0, 1]. More precisely, there exists Xq € (0, 1) such 
that u'(x) < for < x < Xq and u'(x) > for xq < x < 1. 

(ii) u(0) = 0, u(l) = ^. 
( Hi ) For < x < 1 , 

1 , s 1 

_ < u ( x < 

4m 2m 

The proof of Lemma[31is given after the proof of Theoremnbelow. We also need 
the following elementary result from Fourier analysis. 

Lemma 4. For q > 3, < x < tt/2,. 

Px 

< W(x) <— + - < 1.218. 

7T 3 

Proof. As the factor sins in W{x) — ^ J Q s ? g 3 " as is increasing, a standard ar- 
gument in the analysis of the Gibbs phenomenon shows that for < x < ir/2, 
< W(x) < 1 Cj q ds = ^ C si " f 7 , dt. But for < t < tt/2, g i-> gsin(t/g) 

— \ ' — 7T JO sins 7T JO gsin(t/q) — — ' ' ^ * V / 1/ 

is increasing, and so for g > 3 and < a; < 7r/2, < P^(a;) < -| Jj^ 3 si s n " t / 3 ) di = 
^ + |. □ 

Assuming Lemma[31 we now prove Theorem^ By (11711 . integrating by parts 
and using Lemma [IJii), 

r 7r / 2 



4 risings 

y(g) = — / u(coss)as 

7T Jo sin s 



/■tt/2 

2(W(s) m(coss))|q /2 + 2 / W(s)u'(coss) sin sos 
Jo 

/■arccosa^o r 77 /^ 

2( + )W(s) u' (cos s) sins ds. 

«/ «/ arccos a:n 
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Thus, by Lemma |2| and Lemma 

f> arccos xq 



/•arccos xo 

Q(q)<2 / W (s)u' (cos s) sin sds 

Jo 

< 2 ■ 1.218- (u(l) - u(x )) 

< 2 • 1.218 



f\ 1 \ _ 3 1.218 
V 2m 4m /2m 
On the other hand 

,7T/2 

Q(q) > 2 / IL^s) u' (cos s) sins g?s 

•/ arccos xo 

1 1918 

> 2 • 1.218 • (u(x ) - u(0)) > — 

v ' 2 m 

and thus 

3 1.218 1.827 

Q(q) < x = ■ 

2 m m 

Recalling the definitions of <5' m_1 ' and i?!™ -1 !, we have for 1 < i, j < m — 1 



\(Q [m - 1] B [m -%,\ < \Y / Q{n-2i + 2l)2m[,_ 1 ) < 2- 1.827 ■ ^ 
and hence 



J ~ I 



KqIt"-!]^-!]^! < 2 • 1.827 • 1 < i, j < m- 1, 

where L is the rank 1 matrix with entries Lij = 5Zi=o (")' independent of i. Hence 
L has only 1 non-zero eigenvalue Ai(L) and we find 

m— 1 m— 1 



r(L) = A 1 (L) = X;^=EE( I _ 1 



fc— 1 A;— 1 / — 1 

^2m — 1\ m /2m- 1\ o m /2m - 1 

I J l O z 77? — o ^ ( 



1=0 



E, , _ \ / «/t — i \ m / ^m — i \ „■ 



Z 7 2 \ m-l / ~ 2 V m- 1 



In the second last step, we have used the elementary formula preceding (6.7) in 

innu. 

Assembling the above results and recalling the definition of K^ m x l, we obtain 
for 1 < i,j < m — 1, 

= ^fv \(Q lm - 1] B [m - 1] h\ < 

i y i m(2mj! 1 1 J 

where 

(m!) 2 
m(2m) 

and by JTSJ), the only non-zero eigenvalue of K' satisfies 

A 1 (X0=r(^) = 2-1.827-^ T r(L) 

m(2m)! 

„ , nnw (™D 2 m/2m-l\ 1.827 , . , 

< 2- 1.827 y ' s , — , = — ^- = 0.9135 < 1. 

m(2m)! 2\m-l 



Kij = 2 • 1.827 • ^- L« , 1 < i, j < ro - 1, 
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Thus by Lemma |3 

det(l - K [m - 1] ) > det(l - K') = 1 - \ X {K') > 0.0865. 
This completes the proof of Theorem ^ 

Remark 2. Using Lemma [21 the calculations in |DG1| also yield a quantitative 
version of Q but with a weaker bound. As above, we estimate ji™ -1 ! elementwise 
with a rank one matrix so that we can estimate the determinant by estimating the 
only nonzero eigenvalue. We note that we cannot use jDGll (6.22)] (the matrix in 
(6.22) is not rank one). For "small" m we use jDGll (6.15), (6.16)], and for "large" 
m we use fTxTTI (6.55), (6.21)]. We claim that 

(20) detT [m - 1] > 0.02, m > 2. 

This estimate is not optimal, but we could not strengthen it compared to © by 
the methods in jDGl| . To prove (J20JI for 2 < m < 46, we note that the RHS in 
|DG1I (6.16)] is < 0.98 for m in this range. (Note that our Q(q) and I(q) in [DTTT] 
are related by I(q) = mQ{q) - 1 and hence \Q(q)\ = £|1 + I(q)\ < £(1 + |/(?)|).) 
To prove (|20|l for m > 47, we set S = 0.04 and consider q = 3, 5, • • • 4m — 5 in the 
regions ^ ^/ m + 1 / 2 < ^ _ ^ anc [ 1 ^/m+i/2 > ^ ^ se p ara tely. In the former g-region, 
by |DTTT1 (6.21)], |1 + I(q)\ < 1 + \I(q)\ < 1.96. In the latter g-region, substituting 
< -jzr$ m |L)G1I (6.55)], we note that the resulting estimate on |1 + I(q)\ 



multiplied by (± - J^)! 2 2m ~ 2 ), is < 0.98 in fact for m > 44. These facts together 
with Lemma prove (cf. pTTI (6.56), (6.57)]). 

It remains to prove Lemma El A straightforward computation using (|15|) (T^ and 
(|16H shows that u is a solution of the equation 

2 

(21) x(l - x 2 )u' - Amu 2 + (2(m + l)x 2 + 1 - 2m)u - — = 0. 

2m 

Moreover as /i(x) > 0, u is smooth. By H15|)fiiL and by differentiating 1|21|) . we 
find, 

u(0) = 0, u'(0) = 0, u"(0) ' 



/ 99 n m(2m - 3) 

[ ' 1 2 1 

u(l) = , u'(l) = - + . 

V ' 2m' W 3 3m 

Now observe that at a point < x < 1 where u'(x) = 0, we cannot have 4m(m + 

l)u(x) — 1 = 0, i.e. u(x) — 4m (^ +1 ) ■ Indeed, substituting these values into i|2TJl. we 

find —1 + (1 — 2m)(m + 1) = 0, which is a contradiction. Next we show that 

(23) 

, .. . s ii/ \ (Amim + l)u(x) — l) 2 

(u'(x) = for some < x < l) =^> u" (as) = V , , ; , .. . 

v w ; y J m(l - 2mu(x))(l - 4mu(x)) 

Indeed, differentiating l|21|) , we find for such a point x 
,24, „»(*) . '-^It')"' 1 ' - 

Setting u'(x) = in 121|) and solving for (1 — x 2 ) in terms of u(x), we obtain 
^ 1 _ x 2 = _(l-2mu(x))(l-4mu(x)) 



4m(m + l)u(x) — 1 
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Note that by the above argument, the denominator in (|25(l is non-zero: also the 
numerator is non-zero as 1 — x 2 ^ 0. Substituting Ij25(l into l|24|l we obtain (|23|l . 
Furthermore, the calculation shows that if u'{x) — for some < x < 1, then 
u"(x) is (finite and) non-zero. 

From (|22|l we see that for small x > 0, u(x) < 0. As u(l) > 0, there must be at 
least one point x G (0, 1) where u{x) = 0. But it follows from l|21(l that if u{x) = 0, 
x G (0, 1), then u'(x) = 2m (i-x 2 ) > ^- Hence u crosses the level zero at a unique 
point Xi G (0, 1). Next suppose that u'{x) = for some x G (0, X\). But then by 
(|23|l . u"(x) > as u(x) < 0. Thus any critical point for u(x) in (0,xi), must be a 
local minimum. As u(x) clearly has a minimum on (0,xi), it follows that it has a 
unique minimum at xq G (0, say, and no other critical points on (0,xi). Thus 
u'(x) < for < x < xq, and u'(x) > for xq < x < x\. 

Next we show that 

(26) < u(x) < -J— for xi < x < 1. 

2m 

Indeed, if u(x) = ^ for < x < 1, then from (EU we find w'(x) = ^L±I > o. But 
we know from l|22[l that it(l) = 5=-, it'(l) > 0. Hence u(x) cannot cross the level 
^ for < x < 1. This proves (|2l)|) . 

To complete the proof that u is unimodal we show that u'(x) > for x\ < 
x < 1. Suppose 1/(2:2) < for some xi < x 2 < 1. Then as u(xi) = and 
^(2:2) < u(l) = 2^7;, there must exist x\ < x 2 < x 2 and x 2 < x 2 < 1 such that u 
has a local maximum at x 2 and a local minimum at x 2 . By Q28JI. we must have 
u(x 2 ) > and u(x 2 ) < This implies, in particular, that u(x) crosses the 
level j^- at at least one point x# g (x 2 ,x 2 ) such that u'(x#) < 0. But by 1)2 HJ. 
u(x) = 7^, < a; < 1, implies u'(x) = ^ > 0, which is a contradiction. Thus 
it'(x) > on (xi, 1). On the other hand if u'(x3) = for some x\ < X3 < 1, then by 
(|23|l . u"(x^) 7^ and so u'(x) changes sign in a neighborhood of X3, contradicting 
u'(x) > on (xi,l). Thus u'(x) > for all x% < x < 1. This completes, in 
particular, the proof of part (i) of Lemma |3I 

It remains to show that u(x) — u(x;m) > for m > 2, x G [0, 1]. It turns 

out that ,x = x m = 1J plays a distinguished role. More precisely, as we now 
show, 

(27) u(x m ) > =*> (Vx) > for all x G [0, 1]Y 

4m V 4m / 

To see this, suppose u(x) = — j^ for some x G (0, 1): then from (|21(l we obtain 

(m + 2)x 2 — (m — 1) 



(28) u'[x) 



2mx(l — x 2 ) 



Suppose u(x m ) > — 7—. If u(x) < —4^- for some < x < x m , then clearly 
u(x # ) = i/(x # ) > for some point x* G [x,x m ). But by (2SJ), u'(x # ) < 0, 

which is a contradiction. Similarly if u(x) < —4^7 for some x m < x < 1, there 
must exist a point x# G (x m ,x] such that u(x#) = —4^, u'(x&) < 0. But this 
contradicts lt2"8|) as above. This proves (|2"7Jl . 
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To complete the proof of LemmaOH we must prove u(x m ) = u(x m ;m) > 
m > 2. Set s = 1 — x. From (|15|l fiii'). we obtain 

, , . 4m(l - s) 2 ™" 1 f (1 - r)" 2m , 4m(l - s) 2m /" s (1 - t)" 2to , 
/i(a;) = \ ' / v ; ; dr < \- '— -= / ^ ^ dr. 

Using the elementary inequality < e T ~ s for < r < s < 1, we find 

, , , 4me~ 2ms /* a e 2mT , 4me-^ 2 ^ A 2 „ 

Alia;) < 7- 77 : — = / — =- dr = o o — e dX 

where 

/i = Vzms = \j2m(\ — x). 
In order to prove u(x m ) > —-r—, m > 2, we see that it is sufficient to show that 



+ 1 „ 



> for /i = /Lt OT = y2m(l - x m ). 



By the inequality (1— #-) 1 > 1 + j^, and the elementary fact that 1 < fj, m < \/3, 
m > 2, we see that it is sufficient to show 

F(n) > — G(n) for 1 < n < \/3 
m 

where 

F(^) = /ie^ 2 + 2(1 - /i 2 ) e x ' 2 dX, G(ji) = ^ (/ze^ - e A ' dx) . 

Jo 2 V J I 

But G(/i) is clearly increasing and so it is enough to show 

(29) F(u) > ^y^l for 1 < fx < V3. 

m 

Differentiating F(/j,) we find 

F(l) = e, F'(l) = 3e-4/ e A ' rfA > 2.304 > 

Jo 

/•1 2 

F"(l) = 2e - 4 / e A2 dA > -0.415 
Jo 

F'"(/i) > for > 1. 
Thus for 1 < < V3 

j? ff (i \ n 41 ^ — 

> F(l) + F'(l)(/i - 1) + — l) 2 > e — — (\/3 — l) 2 > 2.607. 

On the other hand G(V$) < 41.3, and if we choose m so that 2.607 > then (2§J| 
will hold. Clearly m > 16 satisfies this inequality. We conclude that u(x m ) > —j^ 
for m > 16. On the other hand, using Maple (only sums and products are involved), 
we find from (JTJ, 1161) 



min (u(x m ) + —) > 0.0129 > 0. 

2<m<15 V 4m/ 



This completes the proof of Lemma |21 and hence Theorem ^ 
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